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ABSTRACT 

If the dark matter consists of supersymmetric particles, y-ray observatories such as the Large 
Area Telescope aboard the Fermi satellite may detect annihilation radiation from the haloes of 
galaxies and galaxy clusters. Much recent effort has been devoted to searching for this signal 
around the Milky Way's dwarf satellites. Using a new suite of high-resolution simulations of 
galaxy cluster haloes (the Phoenix Project), together with the Aquarius simulations of Milky- 
Way-like galaxy haloes, we show that higher signal-to-noise and equally clean signals are, 
in fact, predicted to come from nearby rich galaxy clusters. Most of the cluster emission is 
produced by small subhaloes with masses less than that of the Sun. The large range of mass 
scales covered by our two sets of simulations allows us to deduce a physically motivated 
extrapolation to these small (and unresolved) masses. Since tidal effects destroy subhaloes in 
the dense inner regions of haloes, most cluster emission is then predicted to come from large 
radii, implying that the nearest and brightest systems should be much more extended than 
Fermi's angular resolution limit. The most promising targets for detection are clusters such 
as Coma and Fornax, but detection algorithms must be tuned to the predicted profile of the 
emission if they are to maximize the chance of finding this weak signal. 
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1 INTRODUCTION 

Annihilation radiation at y-ray frequencies offers one of the most 
exciting prospects for non-gravitational detection of cold dark mat- 
ter, and is ex pected if the dark matter consists of supersymmetric 
particles (e.g .|Berezinskv et al. 1994, 2003; Bergstrom et al. 1998; 
Stoehr et al. 2003.; .Koushiappas et al.. .2004.; .Colafrancesco et al.. 



atoenr et al.l I2UUA; .Kousmappas et al.. .zUU4.; .Lolatrancesco et al.. 
20071; iDiemand et al.l |2007|; JKuhlen et alj |2008|; IPieri et alj 2008 : 



Springel et al. 2008a; Strigari 
Ackerrnann et al. 2010; Zavala 
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Much effort is being 
devoted to searching for this signal around the Milky W ay's dwarf 
comp anions, in particular using the Fermi satellite CAbdo et al.l 
I2OIOI) . 

Predictions for the properties of the annihilation radiation rely 
on a detailed understanding of the structure of cold dark mat- 
ter haloes which can be gained only through high-resolution nu- 
merical simulations of halo formation. The structure of galaxy- 
mass cold dark matter haloes has been inve s tigated in consider- 
able depth (e.g. lOiemand et all I200I I2OO8I ; iKuhlen et al] I2OO8I ; 
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Springel erai]|2008Jlb t I Anderson et alj201(il;lKamionkowski et al.l 
20101) showing that the radial distribution of low-mass subhaloes, 
and thus of annihilation radiation, is much less centrally concen- 
trated than that of the dark matter as a whole. In the Milky Way, 
this results in the dominant subhalo contribution to the annihi- 
lation radiation coming from large galactrocentric distance and 
so appearing almost uniform across the sky to an observer on 
Earth (Springel et al. 2008a). This same effect causes the annihi- 
lation radiation from an external galaxy cluster to appear much less 
centrally concentrated than the distribution of galaxies. As we show 
below, this has significant implications for the optimal strategy for 
detecting the annihilation signal. 

In this paper we present some of the largest high-resolution 
simulations of cluster haloes to date (the Phoenix Project) and use 
them to investigate the detailed structure of the dark matter distri- 
bution in clusters and its halo-to-halo variation. We use these data, 
together with data from the Aquarius set of galaxy halo simulations 
JSpringel et alJuOOSbl) . to predict the expected y-ray annihilation 
radiation from cluster haloes which we compare to the expected 
annihilation radiation from giant and satellite galaxy haloes. 

As we were completing this work, IPinzke et alj ( 1201 ih 
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and ISanchez-Conde et alj j201lh posted preprints investigating, 
amongst other things, the y-ray annihilation radiation expected 
from galaxy clusters. The luminosity and spatial distribution of 
this radiation depend sensitively on the properties of surviving 
dark matter subhaloes down to the limiting mass of the cold 
dark matter power spect rum, which may be in the r ange 10^* 
to 10~'^ M^ l^ iofmann er^l200ll : iGreen et alj|200^ . For their 
analysis, Pinzke et alj ( 1201 ih relied on an extrapolation of scal- 
ings based on published results for simulations of galactic dark 
matter haloes, including th ose of the Aquarius Project, while 
Sanchez-Conde et alj (1201 ih extended the semi-analytic model of 
Kamionkowski et alj ( l201(]f) . rescaling relevant model parameters. 
Combining the Phoenix and Aquarius si mulations we test ex plic- 
itly the validity of the scal ings used by iPinzke et al.l ( 1201 ih and 
ISanchez-Conde et al.l ( 1201 If) , we investigate their underlying phys- 
ical basis, and we thus construct a more robust (though still uncer- 
tain) framework for extrapolation. For the most part, our results are 
in agreement with those of [Pinzke et alj ( 1201 ih , but not with those 
of lSanchez-Conde et al.l ( 1201 lb who infer a much weaker contribu- 
tion from subhalos to th e total annihilation radiation from clusters 
than IPinzke et alj ( 1201 ih or us find. In this study, we also present 
an estimate of the expected signal-to-noise of the annihilation radi- 
ation from nearby clusters and compare it to that from nearby dwarf 
and giant galaxies. 

The outline of our paper is as follows. Sections 2 and 3 give 
brief descriptions of our simulation suite and of a model for cal- 
culating the annihilation flux and its signal-to-noise in an idealised 
experiment. In Section 4, we discuss our results and their implica- 
tions for dark matter detection. 



2 NUMERICAL SIMULATIONS 

The new dark matter simulations analysed in this study come from 
the Phoenix Project (Gao et al. 2011, in preparation). We supple- 
ment them with previous high-resolution simulations of galactic 
halos from the Aquarius Project carried out by the Virgo Con- 
sortium (Springel et al. 2008a,b). Starting from initial conditions 
appropriate to the ACDM cosmology, both sets of simulations in- 
tegrate the orbits of large numbers of particles using the Gadget- 

3 N-body code (see Springel et al. 2008a). The cosmological 
parameters adopted for both the Aquarius and Phoenix project s 
are those of Virgo's Millennium Simulation dSpringel et al.ll2005r) : 
n,„ = 0.25, flyv = 0.75, CTg = 0.9, Us = 1, and a Hubble constant 
Hq = 100/jkms^' =73kms^'Mpc^'. These were close to the best 
fit v alues derived from t he first year of data from the WMAP satel- 
lite dSpergel et al.ll20031) but are not consistent with the parameter 
ranges found through analysis of the seven-year WMAP data to- 
gethe r with other large-scale structure observations ( Komatsu et al.l 
l201lh . The small offset is, however, of no consequence for the top- 
ics addressed in this paper. 

For the Phoenix Project, we have carried out a suite of ex- 
tremely high resolution simulations of the dark matter distribu- 
tion in galaxy clusters. This suite consists of nine cluster-size dark 
matter haloes with masses in the range [5 — 20] x 10 /?^'M,t;. 
These were selected at random from the Millennium Simulation 
and resimulated at various numerical resolutions. The largest of 
these "Phoenix" simulations, labelled Ph-A-1, represents the dark 
matter with I.O x 10 particles within r2oo. the radius at which the 
enclosed mean density is 200 times the cosmic critical density. It 
has a particle mass of 6.4 x W^h^ Mq and a Plummer-equivalent 
force softening of 0.15/i^'kpc in comoving coordinates at all times. 



This particular cluster has also been simulated at four lower reso- 
lution levels (producing Ph-A-2 to Ph-A-5) in order to assess how 
resolution affects inferences about cluster structure. At the next- 
to-highest resolution level (~ 1.3 x 10 particles within ^200), we 
have simulated an additional eight clusters (Ph-B-2 to Ph-I-2) with 
a particle mass of about 5 x IO^/j^'Mq and a force softening of 
0.32A^'kpc in order to quantify the cluster-to-cluster variation in 
dark matter properties. We will present details of the Phoenix sim- 
ulation suite in a forthcoming paper (Gao et al. 2011, in prepara- 
tion). 



3 RESULTS 

3.1 The total cluster surface brightness 

The total y-ray annihilation luminosity of a dark matter halo is 
the sum of contributions from the smooth main halo, from re- 
solved subhaloes, and from unresolved subhaloes. (Caustics and 
tidal strea ms make a negligible contrib ution to the annihilation 
luminositv iVogelsberger & White! ( 120111) .) If the density distribu- 
tion in the inner regions of the smooth main halo and the re- 
solved subhaloes_Meassumed_to be adequately fit by the NFW for- 
mula dNavarro et al .11 19961 119971) their emission in tegrals [p dV 
can be estimated simply as 1.23 V'^ax/(G 'max) dSpringel et al.l 
uOOSal) . Here Vmax is the maximum circular velocity of the halo 
or subhalo and rmax the radius at which this maximum circular ve- 
locity is reached. 

In Figure[T] we show the azimuthally averaged surface bright- 
ness profile for Ph-A-1, split into the components due to the smooth 
dark matter distribution and to subhaloes resolved down to four 
mass thresholds differing by factors of ten. The subhalo component 
is clearly much less centrally concentrated to the cluster centre than 
the smooth component and its shape appears independent of mass 
threshold as far as can be judged given the noise introduced by the 
finite number of subhaloes involved. The overall level of subhalo 
emission increases steadily as the threshold decreases. The small- 
est subhaloes resolved in Ph-A-1 have masses -^ 5 x 10^ M0, well 
below the masses expected for the haloes of luminous galaxies but 
far above the lower limit for s ubhaloes in a ACDM universe which 
could be as low as lO^'^M© dHofmann et al.ll200ll : lBertone et al.l 
1 20051) . Considerable extrapolation is thus necessary in order to esti- 
mate the total subhalo emission. Note that even at the Ph-A-1 reso- 
lution threshold of 5 x 10^ M0, the surface brightness is dominated 
by the subhalo component at radii greater than 200 kpc. 

To calibrate the extrapolation to lower subhalo masses we 
combine results from our nine Phoenix simulations with results 
from six higher r esolution simulations of galaxy haloes from the 
Aquarius Project dSpringel et alj2008al lbl). Figure 2 shows the total 
annihilation luminosity per unit halo mass (M200) and per decade 
in subhalo mass from subhaloes with masses ranging over 7 orders 
of magnitude, from 10^ to 10 Mg. In the overlap region between 
10 and 10 M0, the Phoenix and Aquarius results agree to about 
30%. This is within the scatter expected given the finite number of 
realizations (illustrated by the shaded area) and the roll-off as sub- 
halo mass approaches 1% of the parent mass. Well away from these 
cutoffs, the shape of this curve is very similar to that of the halo lu- 
minosity per unit mass expected for the Universe as a whole, shown 
as the dashed magenta curve in Figure 2. This reflects the fact that 
the luminosity is domin ated by subhaloe s in the outer regions which 
were accreted recently dCao et al.ll2004h and so have similar lumi- 
nosities and abundance per unit mass (apart from a small bias cor- 
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Figure 1. Surface brightness profiles from dark matter annihilation for var- 
ious components of the Ph-A- 1 simulation of a rich galaxy cluster. Sur- 
face brightness is given in units of annihilation photons per cm^ per second 
per steradian for fiducial values of lOOGev for Hip, the dark matter parti- 
cle mass, and 3 x lO^^^cm^s^' for (ov), the thermally averaged velocity- 
weighted annihilation cross-section, assuming Ny = l photons per annihila- 
tion. This surface brightness scales as Ny{av)/m^ Projected radius is given 
in units of kpc. The red line shows radiation from the smoothly distributed 
dark matter within the main component of the cluster. The ragged blue dot- 
ted lines show radiation from resolved dark matter subhaloes with masses 
exceeding 5 x lO', 5 x 10**, 5 x lO' and 5 x lO'^M^j (from top to bottom). 
Extrapolating to mass Hmits of 10"* and 10"'^M,;i as discussed in the text 
gives rise to the smooth blue curves. The purple dashed lines show the re- 
sults of summing smooth and subhalo contributions. 



rection of 1.5) as the haloes in a representative volume of the Uni- 
verse. Thus, we can use analyti c predictions for the abundance and 
conce ntration of field haloes ( She th & TormenI 120021 : iNeto et al.l 
l2007i) to extrapolate our simulation results to much lower sub- 
halo masses. The upper blue curves in Figure 1 show the resulting 
predictions for minimum subhalo masses of 10"^ and 10" '^M^, 
respectively. The most uncertain part of this extrapolation is the 
assumption that halo concentration continues to increase towards 
lower masses in the same way as measured over the mass range 
simulated so far. This assumption has not yet tested explicitly, and 
has a very large effect on the results. For example, if all (sub)haloes 
less massive than 10 M,^ are assumed to have similar concentra- 
tion, then the total predicted emission from subhaloes would be 
more than two orders of magnitude below that plotted in Figure 1 
for an assumed cut-off mass of 10~ M©. 

With our adopted concentration scaling, subhaloes dominate 
the surface brightness beyond projected radii of a few kiloparsecs, 
as may be seen in Fig. 1. Surface brightness is almost constant be- 
tween 10 and 300 kpc, dropping by a factor of two only at 460 kpc. 
At the virial radius of the cluster (r2oo = 1936 kpc), the surface 
brightness of the subhalo component is a factor of 14 below its 
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Figure 2. Annihilation luminosity (in arbitrary units) from subhaloes lying 
within r2oo per decade in subhalo mass and per unit halo mass (M200) for 
the Phoenix and Aquarius simulations. The level- 1 simulations are shown 
by the black (Phoenix) and red (Aquarius) lines and the medians of the nine 
Phoenix and six Aquarius level-2 simulations by the thick blue and orange 
lines respectively. The full scatter in each set of simulations is indicated by 
the shaded areas. The dashed magenta fine gives the predicted annihilation 
luminosity density per decade in halo mass from the cosmic population of 
dark matter haloes. 



central value. Within this radius the luminosity from resolved sub- 
haloes in Ph-A-1 is more than twice that from the smooth halo, 
even though these subhaloes account only for 8% of the mass. Ex- 
trapolating to minimum subhalo masses of 10"* and 10"'^ Mq 
the subhalo excess becomes 718 and 16089 respectively. These 
boost factors substantially exceed the equivalent factors predicted 
for the galaxy haloes of the Aquarius Project. This is because of 
the additional high-mass subhaloes which contribute in the cluster 
case (see Figure 2) together with the lower concentration of cluster 
haloes relative to galaxy haloes, which reduces the emission from 
the smooth component. Note, the boost factor for the Aq-A-1 ob- 
tained with the extrapolation we use here is smaller by a factor of 
2.4 than the value quoted in Springel et al. (2008a). 

For the resolved component, there is significant variation 
amongst the nine Phoenix haloes, but the median value of the total 
boost factor (for a cutoff mass of 10" Mp,) is 1 125, which, for the 
reasons just given, is about twelve times the median boost factor we 
obtain by applying the same method to the Aquarius haloes. Com- 
paring these results suggests that the ratio of subhalo to smooth 
main halo luminosity within r2oo (subhalo "boost factor") varies 
with halo mass approximately as 



b(M2iyo)=Lsub/Ln 



1.6xlO"-'(M2oo/M0) 



0.39 



(1) 



The total luminosity of a halo is therefore L^ 



= (1 +fe)i'main. 

In addition, the 

projected luminosity profile of the subhalo component can be well 
approximated by 



•Ssub('') = 



16i(M200)irr 



1 



itln(17) 4)0-^16r2 



(2) 



These formulae will be used to estimate dark matter annihilation lu- 
minosities and surface brightness profiles for haloes with different 
masses in subsequent sections. 
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Figure 3. Left panel: Predicted surface brightness profiles of anniliilation radiation (in units of annihilation photons per cm^ per second per steradian) for a 
dwaii' galaxy (UMall; green line), for the nearest large galaxy (M31; red line) and for a rich galaxy cluster (the Coma cluster; black line). As in Figure 1, 
surface brightness scales as Ny{av)/mj,. Projected radius is given in arc minutes. The inner steeply rising part of each curve is due to smoothly distributed dark 
matter in the main halo, while the shoulder of extended emission is produced by low-mass subhaloes. Each profile is truncated at r2oo, the nominal radius of 
the dark matter halo. Right panel: Estimates of the signal-to-noise ratio within a circular aperture of radius R (in arc minutes). The signal is obtained by direct 
integration of the corresponding curves in the left-hand panel and the noise is obtained as discussed in the text. S/W scales as NyB^''^{av)/m^, where Bis the 
surface brightness of the background, assumed to be uniform. 



3.2 Surface brightness and signal to noise of galaxies and 
clusters 

Putting together results from the Phoenix and Aquarius projects, we 
can assess the relative ease of detection of cluster, galaxy and dwarf 
satellite haloes. In the left panel of Figure 3, we show predicted 
surface brightness profiles for three of the most promising candi- 
dates, the Coma cluster of galaxies, the Andromeda Nebula (M31) 
and the dwarf satellite galaxy, Ursa Major-ll (UMa-lI), assuming a 
minumum subhalo mass of 10^ Mq. We represent Coma and M31 
by scaling Ph-A-1 and Aq-A-1 t o the appropriate virial ma sses, 
M200 = 1.3 X 10'^ Mq f or Com a (Reipric h & B6hringedl2002l) and 
1. 8 X lO'^Mfl forM31 JLi & White„200g ). We model UMa-II as 
in lSpringel et alJfeOOSal) (including the contribution from substruc- 
tures - the 'subsub' component). At projected radii below 2 arcmin, 
M3 1 is about twice as bright as UMa-Il and both are substantially 
brighter than Coma. However, at 20 arcmin the surface brightness 
of Coma exceeds that of M3 1 by a factor of 4 and that of UMa-Il 
by about a factor of 6. Beyond about 70 arcmins, M3 1 is again the 
brightest object. 

For y-ray telescopes like the Fermi Large Area Telescope 
(LAT), the detectability of extended objects depends on their con- 
trast relative to the diffuse background. As a simple indicator of 
signal-to-noise (S/N), in the right panel of Figure 3 we estimate 
the signal within a circular aperture from the enclosed luminosity, 
and the noise as the square root of the background counts, assumed 
to be 5 X A X f, where B is the background count rate per unit area, 
A is the area of the aperture in square arc minutes, and t is the 



exposure time. (This assumes that the background is uniform and 
larger than the signal, which may not be the case for the small- 
est apertures.) For the dwarf galaxy UMa-lI, the effective 5/A' is 
almost independent of aperture for radii less than 10 arcmin, but 
drops dramatically at larger radii. In contrast, the 5/A' for Coma 
rises steeply with increasing aperture to a peak at a radius of about 
30 arcmin, significantly larger than the few arcmin resolution of 
the Fermi-LAT at energies ~ 10 GeV. For M31, the effective S/N 
has a minimum on this scale and has maxima on scales of one and 
300 arcmin. In this simple set-up the maximum achievable 5/A' ra- 
tios for Coma and M3 1 exceed that for UMa-lI by about a factor of 
3. 

In practice, realistic experiments will find it difficult to achieve 
these theoretical 5/A' values for very large apertures. Systematic ef- 
fects due to variable backgrounds and difficulties in masking bright 
sources make background correction significantly easier for small 
apertures. M31 is a particularly difficult case because of its very 
large angular size, its low galactic latitude, and confusion from 
other y-ray sources in its inner regions. The Coma cluster is signif- 
icantly more promising because it lies close to the North Galactic 
Pole and appears 10 times smaller on the sky. On the other hand, an 
overly small aperture, corresponding for example to the few arcmin 
resolution of the Fermi LAT instrument at about lOGeV, would 
miss a large fraction of the signal in Coma and other nearby galaxy 
clusters. For a uniform backgr ound, the optimal filt er has a shape 
similar to the predicted profile JSpringel et alj2008al) shown in Fig- 
ure 3 and represented by equation (2). 
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Object Name 


Half-light radius 


Distance 


M200 


L 


F = L/(4jt£/2) 


S/N 




[arcmin] 


[Mpc] 


[Me] 


L^mwJ 


L^mwJ 


[(S//V)n,w] 


AWM7 


35.5 


67.0 


4.2 X 10'* 


7.1x10" 


3.2 X 10-" 


6.8 X 10-3 


Fornax Cluster 


84.1 


17.5 


1.0 X 10"* 


1.2x10" 


8.0 X 10-" 


7.3x10-3 


M49 


59.6 


18.2 


0.4x10'* 


3.9 X 103 


2.4x10-" 


3.1 x 10-3 


NGC 4636 


52.6 


17.4 


0.24x10''* 


2.1 X 10^ 


1.4x10-" 


2.0x10-3 


Centaurus (A3526) 


40.1 


50.5 


2.6x10"' 


3.9x10" 


3.1 X 10-" 


5.8 X 10-3 


Coma 


36.1 


95.8 


1.3 x lO'"' 


2.9 X 10^ 


6.4 X 10-" 


1.3 X 10-2 


Draco 


16.4 


0.082 


N/A 


5.2 X 10"^ 


1.6x10-5 


6.3 X 10-" 


UMal 


18.4 


0.066 


N/A 


4.3x10-3 


2.0 X 10-5 


7.5 X 10-" 


Leol 


4.4 


0.25 


N/A 


3.5 X 10^3 


1.2 X 10-* 


8.2x10-5 


Fornax dwarf 


5.9 


0.138 


N/A 


2.0x10-3 


2.2x10-6 


1.5x10-" 


LeoII 


2.5 


0.205 


N/A 


8.5 X 10-* 


4.1 X 10-' 


3.1x10-5 


Carina 


4.6 


0.101 


N/A 


7.1 X 10-" 


1.4 X 10-6 


1.0 X 10-" 


Sculpt 


13.2 


0.079 


N/A 


3.2 X 10-3 


1.0 X 10-5 


4.9 X 10-" 


Sext 


3.3 


0.086 


N/A 


3.0x10-" 


8.3 X 10-' 


6.1x10-5 


UMall 


28.8 


0.032 


N/A 


2.6x10-3 


5.2x10-5 


1.3x10-3 


Comber 


15.9 


0.044 


N/A 


1.6x10-" 


1.7 X 10-5 


6.8 X 10-" 


Will 


17.7 


0.066 


N/A 


3.9x10-3 


1.8 X 10-5 


7.0 X 10-" 


LMC 


82.5 


0.049 


N/A 


3.8x10-2 


3.3 X 10-" 


3.1x10-3 


SMC 


45.5 


0.061 


N/A 


1.9x10-2 


1.1x10-" 


1.8x10-3 


M31 


351.5 


0.807 


1.8 X 10'2 


1.3 X 10^ 


4.2 X 10-3 


9.3 X 10-3 



Table 1. Principal properties of nearby galaxy clusters, prominent satellites of the Milky Way, and the Andromeda Nebula, M31. The annihilation luminosity, 
L, is given in units of the luminosity from the smooth component of the main Aq-A halo, which we use as a proxy for the Milky Way. The observed flux, 
F, is expressed relative to the flux received by an observer placed 8kpc from the centre of Aq-A. Similarly, the predicted S/N for an optimal filter placed on 
each object is normahzed to the signal-to-noise predicted for a similar filter tuned to the diffuse emission of Aq-A seen from this observer location. For the 
signal-to-noise calculations, we use the optimal filter of Springel et al. (2008a) assuming the background to be the same everywhere and to dominate the signal 
in all objects. 



4 DISCUSSION AND CONCLUSION 

In Table [U we summarise properties of some nearby astronomical 
objects which are relevant for the detectability of their dark mat- 
ter annihilation signal. We consider six galaxy clusters w hich were 
alread y analyzed by the Fermi collaboration JAckermann et al.l 
UOIQ) . thirteen of the known dwarf satellites of our Galaxy, and 
the nearest giant galaxy, M3 1 . 

For the galaxy clusters, distances were taken from the 
NASA/IPAC Extragalactic Databas^j, and virial masse s, M200 
(based on X-ray data), from lReiprich & Bohringeil ([20021). Values 
for Vmax and r,nax were deriv ed assuming an NFW density pro- 
file (Navarro et al. 1 996I . 1 19970 and the mass-concentration relation 
of Neto et al. (2007). We have verified that this relation is consis- 
tent with our simulation data down to the resolution limit of Aq-A- 
1, which is about 10^ Mg. 

Data for dwarf satel lites were taken from the mass models 
of IPenarrubia et alj (|2008!). Their y-ray luminosities are estimated 
from an emission integral based on t he NFW formu l a, f fidV = 
1.23V^3^/(G2rmax)- As discussed in ISpringel et~ZI l2008a) . the 
annihilation signal due to substructures within Milky Way dwarfs 
(the 'subsub' component) is less than that due to the smooth com- 
ponent of their haloes in almost all cases, so we do not consider it 
here. The distance of M3 1 was also taken from the NASA/IPAC 
Extragalactic Database. We base structural parameters for the 
M31 halo on the Aq-A -1 simulation which has a very similar 
mass dLi & Whitj20o3) . 
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We estimate a "best case " signal-to-noise for ea ch object using 
the optimal filter discussed by Springe l et al.l l l200 8a) and assuming 
a uniform background accross the whole sky which dominates over 
the signal in all objects. In this case, the optimal filter has the same 
shape as the signal, and the signal-to-noise can be written in the 
generic form 



S/N = /shape (9/1 /6p,sf) 



tA^fi 



B 



1/2 



F 



(e|+e?sf)'/^' 



(3) 



where F = L/{4Ttd ) is the photon flux, 9/, the ha lf-light radius, 9p sf 
(~ 10 arcmin for Fermi at the relevant energies ( Mic helsonluOOTI) ) 
describes the point spread function of the instrument, / is the inte- 
gration time, Agff is the effective collecting area of the telescope, 
and B is the background count rate per unit solid angle. The func- 
tion /shape i^) encodes the detaile d shape of the emission profile of 



the signal (ISpringel et al.ll2008ah : it is of order unity and depends 
only weakly on the ratio x = 6/,/6psf. 

Using the techniques discussed above, we can compare the 
apparent y-ray luminosities and achievable S/N ratios for galaxy 
clusters with those estimated by Springel et al. (2008a) for dwarf 
satellites of the Milky Way. Results are shown in the Table. We 
find that the brightest nearby cluster, Fornax, is predicted to ap- 
pear 15 times more luminous than the brightest dwarf spheroidal, 
UMall, and 40-50 times more luminous than UMal, Draco or the 
ultrafaint satellite, Wilman-1. However, the Fornax cluster is quite 
extended on the sky, and, as a result, when optimal filters are used, 
the slightly fainter but more compact Coma cluster has a predicted 
S/N ratio 1.8 times larger and ten times that of the most easily de- 
tectable dwarf spheroidal, UMall. Although the Andromeda Neb- 
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ula is predicted to have comparable 5/A', it is not a promising tar- 
get because of the difficulty in correcting for foreground and other 
sources of emission. Note that the 5/A' predicted for both objects 
is still very small compared to that of the main component of the 
Milky Way's smooth halo. Here also, of course, the main problem 
is in separating annihilation radiation from other y-ray signals. 

The Coma cluster thus offers an order of magnitude better op- 
portunity than any Milky Way satellite for detecting dark matter or 
placing limits on its annihilation cross-section. As we have shown, 
for a high resolution experiment like Fermi, the sensitivity for de- 
tecting such radiation will be enhanced by use of a filter which is 
properly matched to the expected extent of the object. For example, 
for the optimal filter, the 5/^ expected for Coma is about 1.5 times 
higher than the 5/A' for a filter based on the point-spread function 
of the Fermi LAT, assuming 10 arcmin for the latter at the relevant 
energies. Detecting annihilation radiation from the Coma or Fornax 
clusters or the placing of robust and stringent upper limits will also 
require careful subtraction of astrophysical sources and an accurate 
estimate of the background. 
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